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ABSTRACT 

We present a measurement of the mean density profile of Ca II gas around galaxies out to ~ 200 kpc, 
traced by Fraunhofer's H & K absorption lines. The measurement is based on cross-correlating the 
positions of about one million foreground galaxies at z ~ 0.1 and the flux decrements induced in the 
spectra of about 10 5 background quasars from the Sloan Digital Sky Survey. This technique allows 
us to trace the total amount of Ca II absorption induced by the circumgalactic medium, including 
absorbers too weak to be detected in individual spectra. We can statistically measure Ca II rest 
equivalent widths down to several mA, corresponding to column densities of about 5 x 10 10 cm~ 2 . We 
find that the Ca II column density distribution follows iVcaii ~ f~ 1A and the mean Ca II mass in 

the halo within 200 kpc is ~ 5 x 10 3 M©, averaged over the foreground galaxy sample with median 
mass ~ 10 10 ' 3 M Q . This is about an order-of-magnitude larger than the Ca II mass in the interstellar 
medium of the Milky Way, suggesting more than 90% of Ca II in the Universe is in the circum- and 
inter-galactic environments. Our measurements indicate that the amount of Ca II in halos is larger for 
galaxies with higher stellar mass and higher star formation rate. For edge-on galaxies we find Ca II 
to be more concentrated along the minor axis, i.e. in the polar direction. This suggests that bipolar 
outflows induced by star formation must have played a significant role in producing Ca II in galaxy 
halos. 

Subject headings: quasars: absorption lines - galaxies: evolution - galaxies: halos - intergalactic 
medium 



1. INTRODUCTION 

The H & K absorption and emission lines have had an 
important role in the early developments of spectroscopy 
and, for about two centuries, have been key to under- 
standing a host of astronomical phenomena ranging from 
the Solar atmosphere to the interstellar medium of dis- 
tant galaxies. 

The labelling convention for the major absorption lines 
traces its origin back to the work of the German opti- 
cian and glass manufacturer Joseph Fraunhofer who, in 
1814, catalogued over 500 lines in the solar spectrum 
and designated the principal dark fe atures with letters 
from A through H (jFraunhoferl 118141 ) . He used the let- 
ter I to denote the blue end of the visible spectrum and 
labelled a selected set of weaker lines with lower-case let- 
ters. Despite being one of the strongest features in the 
solar spectrum, the bluer line of the doublet was never 
labelled by Fraunhofer. His main goal was to measure 
the refractive indices of various materials and was prob- 
ably not interested in using pairs of lines closely sepa- 
rated in wavelength. The nineteenth-century literature 
re fers to this second line as H2, H' or HII. According 
to Pasachoff fe Suerl (J2010I ). the origin of the modern la- 
belling goes back to John William Draper who in 1843 
named the lines as H and k (lower case) in order of 
decreasing wavelength (Draper 1843). The first person 
to use the capital K to designate the second H line is 
Elcuthere Mascart from the Ecole Superieure de Paris 
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in 1864 (|MascartJ fl864). The naming convention was 
straightened out at the first International Astronomical 
Union meeting, held in Rome in 1922. The Committee 
accepted that only some of the mostly century-old sym- 
bols should be retained. It was decided that only A (O2), 
a (0 2 ), D (Na), b (Mg), G (Fe & Ca), H and K were to 
be preserved and the notation Ha and H/3, etc. sh ould 
be adopted for the hydrogen lines <|Hearnshawl [l986) . 

The H & K lines have played a major role in astro- 
physics because of their strength and their location in 
the visible part of the spectrum. They correspond to 
the fine structure splitting of the singly ionized calcium 
excited states Ca II (Ca+) at A = 3934.78 A (K) & 
3969.59 A (H) in vacuum. Being an abundant element, 
log(Ca/H) Q + 12 ~ 6.34 (|Asplund et al.ll2009D . calcium 
is found in a range of astrophysical environments. Cal- 
cium is a refractory element and has ionization potentials 
of 6.11 and 11 .87 eV, for Ca I and Ca II, respectively 
(|Mortonl [2003) . In high-density environments such as 
the disk of our Galaxy, most of the calcium is found in 
dust grains. In low-density environments, it tends to be 
found in Ca III and other higher-ionization states. 

The H & K lines are the only resonance lines in cool- 
star spectra of the dominant ionization stage of an abun- 
dant element that are accessible to ground based obser- 
vations. In 1814, Fraunhofer himself detected these ab- 
sorption lines in the spectra o f various stars. About a 
century later. iHartmannl (|1904[ ) detected H & K absorp- 
tion due to interstellar matter, opening a new window 
on the mapping of the gas distribution in the Galaxy. 
O n extragalactic scales, Ca II a bsorption was detected 
by iBoksenberg fe Sargent! (|l978) in the spectrum of the 
quasar 3C 232. The first detectio ns of the H fe K absorp- 
tion due to galaxy halos followed f Bokscn berg fe SargentJ 
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Fig. 1. — Individ ual Ca II absorber s with known associated galax- 
ies, compiled from York ct al. (2012, green triangles), Richtcr ct al. 
(2011, magenta squares), and Bowcn ct al. (1991, black diamonds). 
Only those beyond impact parameter r p = 3 kpc are shown, among 
which two have r p larger than 20 kpc. 



[1971 IBoksenberg et al.lll98(l Blades etlIIH98lh . Fig- 
ure Q] shows a compilation of the detections reported in 
the literature. In such low-density environments, the 
mapping of Ca II has been limited to the immediate 
vicinity of galaxies, i.e. within about 20 kpc due to the 
weakness of these absorption lines. The first detection 
of calcium H & K absorption by a high-velocity cloud 
(HVC) was reported by West et al. (1985), showing the 
Milky Way halo is enriched with metals. In this paper 
we extend the use of the H & K lines to the circum- and 
inter-galactic context. 

We present a measurement of the mean density pro- 
file of Ca II in galaxy halos out to ~ 200 kpc. We 
are also able to explore its dependence on stellar mass, 
star formation rate (SFR), and azimuthal angle, and 
place strong constraints on both the total amount of 
Ca II in galaxy halos and its origins. Our statistical 
analysis is based on cross-correlating the positions of 
about one million foreground galaxies at z ~ 0.1 and 
the flux decrements induced in the spectra of about 10 5 
background quasars from the Sloan Digital Sky Survey 
(SDSS, lYorket al.ll2000D . Metals in galaxy halos can 
have different origins - enriched outflows driven by star 
formation or active galactic nuclei, gas stripped out of 
satellite galaxies by ram pressure or tidal force or out- 
flows, gas in satellite galaxies, and enriched infalling in- 
tergalactic gas. Probing the distribution of Ca II may 
shed new light on these processes. 

The paper is organized as follows: In Section 2, we 
briefly describe the dataset and the method. In Section 
3, we present the results of the statistical analysis, esti- 
mate the total Ca II, calcium and gas mass, and explore 
their dependence on galaxy properties. Section 4 sum- 
marizes our results. We assume the ACDM cosmogony, 
with Q A = 0.7, Q m = 0.3, and H = 70 km s _1 Mpc -1 . 

2. DATA ANALYSIS 

Our goal is to statistically measure the H & K absorp- 
tion lines imprinted in the spectra of distant quasars by 
the circumgalactic medium (CGM) of low-redshift galax- 
ies and estimate the mean column density of Ca II within 
their virial radius and beyond. To do so we measure the 
mean rest equivalent width, (Wq ), of the Ca II K line 



imprinted in quasar spectra, as a function of impact pa- 
rameter (projected galactocentric distance, r p ) from low- 
redshift galaxies. A similar technique has been used at 
higher redshi f ts to p robe various UV absorption lines by 
ISteidel et al.] (|2010|) with the Ke ck Baryonic Structure 
Survey and IBordoloi et~aT1 (|2011[ ) using the zCOSMOS 
survey (|Lillv et al.ll2007fl . 

When Ca II absorption is in the linear regime of the 
curve of growth, the rest equivalent width is related to 
Ca II column density with 



A^Caii = 1.13 X 10 20 cm" 2 



^ (1) 

/A 2 ' l ' 

in unit of A and the 

Ca II K line is 0.648 

This relation is valid 

Ca II K absorption is 



where both W^ and A are 

oscillator strength / for the 

(|Safronova fc Safronoval 1201 If) 

when the optical depth for 

smaller than unity, which corresponds to -/Vcaii ^ 

10 125 cm~ 2 (assuming the Doppler broadening factor 

b ~ 10kms _1 j3- As shown below, this is the regime 
of interest in the present study. 

Our analysis makes use of data from SDSS DR7 
([Abazaiian et al.lr2009l). We use the value-added MPA- 
JHU catalog ( Kauffmann et ail [2003; B rinchmann et ail 
2004J3, which includes stellar mass and SFR measure- 
ments for 819, 437 unique galaxies with a mean redshift 
z ~ 0.1. We also supplement the catalog with position 
angle and minor-to-major axis ratio (b/a), retrieved from 
SDSS skyserveiQ- To ensure that the induced Ca II ab- 
sorption is not blended with that of our Galaxy, we se- 
lect galaxies at redshifts z ra i > 0.02. We selec t quasars 
as background sources ([Schneider et alJ | 2010D and use 
the improved redshift estimates given bv lHewett fc Wild! 
(20100 The sample includes 107, 194 quasars with red- 
shifts spanning 0.1 < z < 6.5. 

Absorption measurements are based on continuum- 
normalized fluxes: 



R(X) 



F(X) 



t(A) 



(2) 



For the extraction of weak absorption information from 
the noise-dominated flux residuals, accurate estimation 
of the source flux continuum is crucial. We make use 
of the SDSS quasa r conti nuum estimates F cont (X) given 
by iZhu fc Menard (|2012f) . In a nutshell, this analysis 
employed the robust dimensionality-reduct ion technique 
nonnegative matrix factorizat ion (NMF, iLee fc Seungj 
1999: Bla nton fc Roweisl I2007T ) to construct a basis set 
of nonnegative quasar eigenspectra, and fit each ob- 
served quasar spectrum with a nonnegative linear com- 
bination of these eigenspectra. Large-scale residuals not 
accounted for by the NMF basis set were removed with 

4 From high-resolution spectroscopy of individual Ca II ab- 
sorbers, the median Doppler broadening factor b ~ 6km s _1 (e.g., 
IRichter ct al. 2011). Usually there are several unsaturated compo- 
nents induced by the CGM of a single galaxy, and we cannot resolve 
these different components here because SDSS has an instrumental 
resolution ~ 69 km s — 1 . The apparent Doppler broadening factor 
therefore is dominated by the separation of different components, 
which is 3> 10 km s . 

5 http : //www . mpa-gar chlng . mpg . de/SDSS/DR7/ 

6 http : //cas . sdss . org/astrodr7/en/ 

7 http : //das . sdss . org/va/Hewett_Wild_dr7qso_newz/ 
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Fig. 2. — Stacked continuum-normalized spectra of background quasars. In the left panel, the gray lines show the geometric means of 
original residual spectra, while the blue lines show the geometric means of residual spectra convolved with a double Gaussian profile. In 
the right panel, we compare the Ca II signals (dark blue vertical line) with the distributions of the flux residuals (light blue histogram). 



appropriate median filters. For details regarding this 
continuum estimation method we refer the reader to 
iZhu fc Menardl (|2012f ). The continuum estimation re- 
quires robust NMF basis, which restricts the quasar sam- 
ple to those at z < 4.7. In total we are able to use 
84, 533 quasar spectra. Finally, when selecting back- 
ground quasars, we require z qua sar — z m \ > 0.1. 

As we are intersted in estimating absorption equivalent 
widths in the rest-frame of selected galaxies at various 
rcdshifts, we will shift every flux residual R(X) to the 
rest frame of the absorbing material throughout the pa- 
per. To optimally extract information on absorption lines 
we us e a matched-filter technique (e.g., iSchneider et al.l 
1993) to convolve the residuals R(X) with the expected 
absorption line profile P (A): 



equivalent width is given by 



< = -[l-i? c (A K )] 



(5) 



RcW = 1 - 



/[l-.R(A')]P(A'-A)dA' 



/P 2 (A'-A)dA' 
where P(A) is normalized with 



J P(X)dX= 1 



(3) 



(4) 



When measuring the mean Ca II rest equivalent width 
induced by a sample of foreground galaxies, we select 
quasars for which the expected Ca II absorption fea- 
ture does not fall in the Lya forest. We also exclude 
quasars for which the expected Ca II absorption over- 
laps the quasar's C IV, C III], and Mg II emission- line 
regions, although including these systems in the analysis 
has little effect on our results. 

3. RESULTS 

3.1. The average Ca II absorption radial profile 

Our estimate of the mean Ca II absorption is given by 
the cross-correlation between the position of foreground 
galaxies denoted by 5 g (r) and the absorption rest equiv- 
alent width expected in the spectra of quasars, as a func- 
tion of projected separation r p : 



«>M = (S s . <M) 



(6) 



When focusing on a single absorption line, P(X) is chosen 
to be a Gaussian kernel, with a fixed velocity dispersion 
of 3 pixels (~ 207 km s~ ). We note the exact value of the 
velocity dispersion adopted (2, 3, or 4 pixels) has little ef- 
fect on the results. To extract information from the Ca II 
doublet, P(X) is a double-Gaussian kernel with a 34.8 A 
separation and a doublet ratio of 2 and for which A is the 
central wavelength of the blue K line. The convolved flux 
residuals can be used to obtain a non-parametric mea- 
surement of the Ca II equivalent width. The K-line rest 



where the ensemble average is taken at the position of 
every galaxy and includes all background quasars in a 
radial bin centered on r p . Absorption being a multi- 
plicative effect we use a geometric mean for the ensemble 
average, which gives us an arithmetic mean of the cor- 
responding optical depth. However, using an arithmetic 
mean yields similar results, as expected when measuring 
weak absorption lines. Our estimator is inverse-variance 
weighted, using the wavelength-dependent noise given by 
the SDSS pipeline. 
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Fig. 3. — The rest equivalent width of Ca II as a function of impact parameter. The blue solid circles show the measurements using 
spectra convolved with the double Gaussian kernel. The gray open diamonds are measurements using K line only (with a single Gaussian 
kernel). The solid line shows the best-fit power law to the measurements at r p < 200 kpc, while the dotted line is the extrapolation of the 
best-fit profile. The inset shows the measurements and the extrapolation of the best-fit profile on large scales, but in linear scale in j/-axis. 
The vertical dashed line shows the mean 90% light radius of the foreground galaxy sample ((rgo} ~ 8 kpc). 



TABLE 1 
The average Ca II absorption radial 

PROFILE 



r p bin 


Median r p 


■* 'pairs 


(W/ K (CaII)) 


[kpc] 


[kpc] 




[mA] 


(3, 10] 


7 


17 


435 ± 68 


(10,20] 


15 


86 


132 ± 33 


(20, 30] 


26 


131 


59 ±27 


(30, 45] 


39 


410 


40 ±15 


(45, 68] 


58 


1106 


19.5 ±8.7 


(68, 101] 


87 


2991 


16.5 ±5.7 


(101,152] 


130 


7456 


5.2 ±3.4 


(152,228] 


194 


17348 


5.8 ±2.4 



In Figure [2] we show intermediate measurements lead- 
ing to our estimate of (W^)(r p ). The left panel 
shows, for four radial bins, the stacked galaxy rest-frame 
& continuum-normalized quasar fluxes (-R(A)) between 
3800 and 4100 A (gray lines). The vertical lines show 
the expected locations of the H & K lines. The blue 
lines show the convolved residuals (i? c (A)). To esti- 
mate the uncertainty, we measure the distribution of the 
convolved residuals (i? c (A)) over the wavelength range 
between 3800 and 4100 A and estimate its dispersion. 
This is illustrated in the right panel of the figure, where 
the vertical lines indicate the corresponding values of 
<i?c(A K )>. 

We present the radial dependence of Ca II absorp- 
tion around galaxies in Table Q] and Figure [3] The blue 
solid circles represent the rest equivalent width estimates 
(Wq) derived from measurements of the H & K lines 
(-Rc(Ak))- We also show the estimates based on the resid- 
uals of the K line only ((J?(Ak)), with a single Gaussian 



kernel) with the gray open diamonds. The conversion to 
column density is shown in the right y-axis of the fig- 
ure. It is important to note that even if some overlap 
might exist between the angular apertures used around 
different galaxies, the measurements on different scales 
are largely uncorrelated as different galaxies tend to be 
at different redshifts. The inset of the figure shows, using 
a linear scale, the measurements on scales greater than 
200 kpc. Our results show that, on average, Ca II gas is 
present in the halo of galaxies up to scales comparable to 
their virial radius. For reference, the vertical dashed line 
shows the mean value of rgo, the radius encompassing 
90% of the light of a galaxy, which we take as a fiducial 
limit between the disk and halo environments. For our 
sample (rgo) ~ 8 kpc. The median stellar mass of the 
foreground galaxy sample is ~ 10 103 Mq. 

To validate the robustness of our measurements, we 
perform two null hypothesis tests: (1) Using Eq. [5] we 
measure the rest equivalent widths at different wave- 
lengths. The left panel of Figure 0] shows the signif- 
icance of the measurements for 12 randomly-selected 
wavelengths between 3800 and 4100 A (gray diamonds). 
For comparison, the Ca II measurements are shown as 
blue solid circles. (2) We measure the expected Ca II rest 
equivalent widths but, instead of using the corresponding 
background quasars lying within some impact parameter 
from the galaxies, we use random quasars with similar 
redshifts. These results are shown in the right panel of 
the Figure. In both cases our test measurements are 
consistent with no signal and show that the detection of 
Ca II absorption above (Table Q] and Figure [3]) is robust 
and not induced by systematic effects. 

3.2. Total mass of Ca II in the halo 
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Fig. 4. — Null hypothesis tests showing the robustness of the Ca II detection. Left panel: significance of rest equivalent width measurements 
at randomly-selected wavelengths. Right panel: significance of rest equivalent width measurements but using random quasars at similar 
redshifts as those in galaxy-quasar pairs. The blue solid circles show the significance of the Ca II absorption measurements. 



We parametrize the radial dependence of the mean 
Ca II absorption around galaxies with a power-law func- 
tional form: 



«)(r p )=A 



100 kpc 



(7) 



where r p is taken as the median impact parameter in each 
bin. Performing a least-squares fit to the measurements 
at r p < 200 kpc, we obtain: 

A= 11 ± 2mA and a = -1.38 ±0.11. (8) 

We show the best-fit power law with the blue line in 
Figure El where we also show its extrapolation to larger 
scales. The slope shows the mean Ca II distribution is 
steeper than an isothermal profile. Converting (Wq ) to 
mean column density (-/Vcaii) using Eq.[T]and integrating 
the best-fit profile between 10 and 200 kpc, we estimate 
the average total Ca II mass in the halo to be 0: 



(M, 



halo \ 
Call/ 



200 kpc 



= 27T 

' 10 kpc 

= (5.0 ±1.0) x 10 3 M C . 



(Ncaii)(r p ) r p dr p 



(9) 



It is interesting to compare the total Ca II mass in 
the halo with the total Ca II mass in galaxies themselves 
such as the Milky Way. To estimate the total Ca II mass 
in the Milky Way disk, we use the average neutral hydro- 
gen H I surface densit y lOM pc~ 2 (~ 10 21 cm~ 2 , e.g., 
iKalberla fc Kerpl 120091 ). and the typical Ca II-to-H I ra- 
tio Nc,f,.u/Nki ~ 10~ 8 - 9 in the interstellar medium (e.g., 
IGudennavar et all I2012T ) . Integrating the surface den- 
sity out to 15 kpc, we obtain M<¥ a Yi Disk ~ 3 x 1O 2 M , 
an ordcr-of-magnitude smaller than the average total 
Ca II mass in the halo, suggesting more than 90% of 

8 We note this cylindrical mass is obtained using projected col- 
umn density. For spherical mass we need to deproject it through 
inverse Abel's integral with uncertain assumptions of gas distribu- 
tion on large scales. The correction factor is ~ 0.6 and the average 
spherical mass is ~ 3 X 10 3 Mq assuming the power-law distribu- 
tion extends to infinity. Throughout the paper we will stick to the 
observable cylindrical mass. 



Ca II in the Universe is distributed in circum- and inter- 
galactic environments. This dramatic difference, how- 
ever, is likely caused by the difference in the dust deple- 
tion levels rather than in the total amount of calcium in 
different environments. From observations of the Milky 
Way halo, the gas density in the CGM is more than an 
order-of-magnitude lower than in the interstellar environ- 
ment, therefore there is relatively les s calcium depleted 
onto dust grains (e.g., Wakk er fc Math is 2000), resulting 
in a Ca II abundance an ordcr-of-magnitude higher. Fi- 
nally, we note that the total Ca II mass in the halo above 
is obtained by averaging over all the foreground galax- 
ies (in the galaxy-quasar pairs), whose median mass is 
~lO la3 M . 

3.3. HI mass traced by Ca II 

We now attempt to estimate the total amount of 
H I gas traced by Ca II. It is known that in 
cool, dense clouds most of Ca is depleted onto grains 
and Ca II can be the dominant ionization state, 
while in warm, lower-density medium, Ca is less 
depleted but the dominant ionization s t ate is usu- 
ally Ca III (e.g.. iRoutlv fc SpitzeH H952T: iWeltv et al.1 
fl99J iWakker fc Mathisl 1200 0: Ben Bckhti et al. 2012). 
Here we make the assumption that the Ca II-to-H I 
in the halo is similar to those observed for high- 
velocity clouds/Magellanic Stream in the Milky Way halo 
{Nce.il/N-g_i ~ 10^ 7 - 2 . IWJtlt^[200l . In this case the to- 
tal amount of H I traced by Ca II in galaxy halos is about 
several x 1O 9 M . For comparison, the total amount of 
H I in the Milky Way halo, combining the Magellanic 
system (Magellanic Stream, the Leading Arm, Magel- 
lanic Bridge, excluding galaxies themselves, 3 x 10 8 M ), 
LMC/SMC interstellar medium (7 x 10 8 M ),andHVCs 
(3_xlO 7 M ), is estimated to be ~ 10 9 M l|Putman et all 
12012ft . The total H I mass in M31's halo, without count- 
ing g as in M33, is ~ severalx 1O 7 M (e.g., iThilker et al.l 
2004). 



3.4. Dependence on galaxy properties 

We now investigate the relation between galaxy prop- 
erties (stellar mass M*, SFR, and azimuthal angle (f>) 
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Fig. 5. — Stellar mass-SFR diagram of foreground galaxies in 
galaxy-quasar pairs. The green line separates star- forming and 
quiescent galaxies. The blue vertical line (M* = 10 Mq) then 
divides star-forming galaxies into high-mass and low-mass star- 
forming subsamples, and the red vertical line (M t = 10 107 Mq) 
divides quiescent galaxies into high-mass and low-mass quiescent 
subsamples. 

and the amount of Ca II in their halos. Considering 
the bimodal nature of the distribution of galaxies in the 
stellar mass and SFR space, we divide the sample of fore- 
ground galaxies using t he following demarcation relation 
(jMoustakas et alJl2013t l: 

log 10 SFR = -0.79 + 0.65 (log 10 M* - 10.0) , (10) 

where SFR is in M0 yr -1 and stellar mass M* is in 
Mq. This demarcation relation is shown as the green 
line in Figure [5] and allows us to separate star-forming 
galaxies from quiescent ones. To divide each population 
into high-mass and low-mass subsamples, we adopt a 
stellar mass separation of M* = 10 101 and 1O 1O ' 7 M 
for star-forming and quiescent galaxies as shown with 
the vertical lines in Figure \5\ This selection takes 
into account the fact that star-forming galaxies are on 
average less massive than quiescent galaxies. We now 
explore how the amount of Ca II in the halo depends on 
the following parameters: 

• Stellar Mass: we select two samples of galaxies us- 
ing the mass cuts defined above, combining star-forming 
and quiescent galaxies. Our measurements show that 
more massive galaxies have more Ca II in their halos. 
To evaluate the difference quantitatively, we fit the Ca II 
density profile with a power-law function with a fixed 
slope a = —1.38 as found in Eq. [5] (the reduction in 
S/N prevents us from placing robust constraints simul- 
taneously on both the amplitude and the slope of the 
measurements). The best-fit amplitudes are 



^high-mass — 15 ± 2 mA 
^t-Iow — mass = I i 2 mA . 



(11) 



For reference, the median stellar masses of high-mass 
and low-mass populations are 10 105 and 1O 9 ' 8 M . 
Thus, within a fixed impact parameter, the dependence 
of the amount of Ca II on stellar mass is shallower than 



Since the halo of higher-mass galaxies is more extended 
than that of lower-mass ones, we also look at the depen- 
dence on stellar mass by measuring the absorption as 
a function of impact parameter in unit of virial radius, 
fpAvir- To estimate the virial radius r v ; r of a galaxy with 
a given stellar mass, we assume a power-law relation be- 
tween virial mass M v - lr and stellar mass M* : 



M* oc M, 



5/3 



(12) 



where the slope value is taken to be approximately the 
average slope determined from the abundance match- 
ing technique for ga laxies selected betwe en M* = 10 9 
and 1O U M (e.g., iBehroozi et"aHl2010l ). At a given 
stellar mass below lO n M , the host dark matter halo 
mass is statistically independent of morphology (thus 
SFR), as shown by gala xy-galaxy lensing studies (e.g., 
Mandclbaum ct al. 2006). The virial radius is then cal- 
culated via: 



47rp c A c (z) 



(13) 



where A c (z) 



(18tt 2 + 82a: - 39a; 2 ) and x = fi m (z) - 1 



(e.g., iBrvan fc Normanlll998l) . At the median mass of 
the z ~ 0.1 foreground galaxy sample M* = 1O 1O3 M , 
M vir ~ 1O 118 M and r vir ~ 200 kpc, which is the outer 
bound we used to calculate the average total mass of 
Ca II in the halo for all galaxies. The virial radius is 
then estimated to be 



200 kpc 



M, 



1Q10-3M 



o 



1/5 



(14) 



The virial radius-scaled profiles are shown in the right 
panel of Figure [51 Fitting the dependences with a power- 
law, keeping the slope a, we find the best-fit amplitudes 
(at r = 0.5 r v ir) are: 



^high-mass — 10 ± 2 mA 



4' 



low— r 



9 ± 2 mA , 



(15) 



consistent with each other. This shows that the amount 
of Ca II in the halo is, on average, more closely related 
to the dark matter mass than the stellar mass. In other 
words, the results indicate that the total Ca II mass 
within virial radius follows roughly: 



rhalo 
'Call 



M, 



3/5 



(16) 



• Star Formation Rate: we now divide the galaxy 
sample into star forming and quiescent galaxies using the 
demarcation represented by the green line in Figure [5] 
(Eq. flO|) . We show the corresponding column density 
profiles of Ca II in the left panel of Figure [7] Our mea- 
surements show that, on average, star forming galaxies 
are surrounded by more Ca II than quiescent systems. 
We fit the Ca II density profile with a power-law func- 
tion fixing the slope a and find the amplitudes to be: 



^•star— forming — 12 ± 2 mA 
^quiescent = 6 ± 2 HlA . 



(17) 
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Fig. 6. — Dependence of Ca II absorption on stellar mass. The left panel shows that the Ca II absorption is stronger around high-mass 
galaxies (green solid circles) than low- mass ones (magenta triangles). The right panel shows the Ca II absorption at fixed virial radius- 
normalized galactocentric distance. The virial radius-scaled Ca II absorption strengths at different stellar masses appear to be consistent 
with each other within errors. The corresponding lines show the best-fit power-law profiles with a fixed slope of a = —1.38. 
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Fig. 7. — Same as Figure \E\ but for SFR dependence. There is more Ca II around star-forming galaxies than quiescent ones, 
dependence is stronger after normalizing impact parameter with virial radius thus removing the stellar mass dependence. 

In the right panel, we show the column density profiles 
normalized by the estimated virial radius of the galaxies. 
This removes the stellar mass dependence and the SFR 
dependence is slightly stronger with best-fit amplitudes: 



The 



TABLE 2 
Average Call mass in the halo 



Sample 



Astar- forming — 13 ± 2 mA 
^quiescent = 5 ± 2 mA . 



(18) 



This result indicates that either star formation (or a 
physical process correlated with it) is responsible for an 
excess of Ca II in galaxy halos or that Ca II is suppressed 
in the halo of quiescent galaxies. 



Median M* 

[io 10 Mg] 



Median r v 
[kpc] 



/ AT Halo \ a 

\ M Ga.lV 

[10* M Q ] 



High-mass star-forming 


3.3 


221 


9.5 ±1.8 


Low-mass star-forming 


0,1 


113 


4.1 ±1.8 


High-mass quiescent 


9.5 


271 


3.8 ±1.8 


Low-mass quiescent 


2.2 


205 


1.4 ±1.8 


ALL 


2.0 


200 


5.0 ±1.0 


a Average total Ca II mass in the halo 


between 0.05 r v 


r and 1 r v i r 



In Figure [5] and Table [21 we summarize the total 
Ca II mass in the halo (measured between 0.05 r v ; r 
and r V i r using Eq. [9]) for the four galaxy subsam- 
ples. In short, we find that there is more Ca II in the 
halo of galaxies with higher stellar mass and higher SFR. 

• Azimuthal Angle: we now constrain the spatial 
distribution of Ca II as a function of the location of the 
gas with respect to the disk of the galaxy. To do so 
we select edge-on galaxies with minor-to-major axis ra- 



tio b/a < 0.55, and for a given galaxy-quasar pair we 
measure the azimuthal angle (f> of the quasar from the 
galaxy's minor axis (determined from the position angle). 
We then divide the galaxy-quasar pairs into two subsam- 
ples: |0| < 45° (along the minor axis) and \cf>\ > 45° 
(along the major axis) and measure the average Ca II 
absorption for each subsample. Figure |H] shows the mea- 
sured Ca II absorption radial profiles. We find that, on 
average, the absorption is stronger along the minor axis 
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Fig. 8. — Dependence of the total amount of Ca II in the halo 
on stellar mass and SFR. The total amount of Ca II is estimated 
by fitting the measured absorption radial profile with a power- 
law form with a fixed slope of a = —1.38 and integrating the 
best-fit profile between 0.05 r v ; r and r v ; r (10 kpc and 200 kpc for 
a 1O 1O - 3 M galaxy). 

(in the polar direction) than along the major axis (par- 
allel to the disk plane). Fixing the slope a, we find the 
best-fit amplitudes are: 



Apolar = 18 ± 3 mA 
Aiisk = 6 ± 3 mA . 



(19) 

This indicates that there is more Ca II perpendicular to 
the disk and indicate that bipolar outflows driven by star 
formation (or a process correlated with it) may be one 
of the major sources for Ca II. We note that no signifi- 
cant excess is detected on large scales. It is likely that 
the azimuthal angle de pendence itself depen ds on the 
galactocentric distance. IBordoloi et al.l (|2011[ ) reported 
a similar effect by measuring Mg II absorption around 
galaxies and was able to detect an azimuthal dependence 
only within ~ 50 kpc from the galaxies. 

4. SUMMARY AND DISCUSSION 

The Ca II H & K absorption lines have been used 
for two centuries to study a host of astrophysical phe- 
nomena. In this paper we have extended the use of 
these lines to the circum- and inter-galactic contexts by 
measuring the galaxy-Ca II gas spatial correlation func- 
tion. Using spectra of about 10 5 quasars from the SDSS 
DR7 dataset we have measured the mean Ca II absorp- 
tion around one million low-redshift galaxies with me- 
dian M* ~ 1O 1O ' 3 M0 on scales ranging from about 10 
to 200 kpc. This statistical technique allows us to probe 
absorption lines down to rest equivalent widths of several 
milli- Angstroms and to be sensitive to the total amount 
of Ca II around galaxies. 

The detection of absorption features relies on the abil- 
ity to accurately estimate the intrinsic spectral energy 
distribution of the background source. To do so we 
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Fig. 9. — Dependence of Ca II absorption on azimuthal angle. 
The cyan solid circles show the Ca II absorption along the minor 
axis, and the orange triangles along the major axis. There is more 
Ca II along the minor axis. The corresponding lines show the best- 
fit power-law profiles with the slope fixed at a = —1.38. 

have used the continuum estimates presented in Zhu & 
Menard (2012), based on a vector decomposition tech- 
nique. The analysis presented in this paper is based on 
the cross-correlation between galaxy positions and the 
flux decrements in the background quasar spectra. It 
has allowed us to obtain the following results: 

• We detect Ca II H & K absorption out to pro- 
jected galactocentric distance r p ~ 200 kpc (about 
virial radius). This is about an order-of- magnitude 
farther out than previous studies based on indi- 
vidual absorber detections. Our sensitivity allows 
us to probe Ca II column densities down to about 
5 x 10 10 cm- 2 . 

• The mean spatial distribution of Ca II gas around 
galaxies follows a power-law column density profile 
with a slope a = — 1.38±0.11, at r p < 200 kpc (see 
Eq. [7|), which is steeper than an isothermal profile. 
We estimate the total Ca II mass in the halo within 
200 kpc to be ~ 5 x 10 3 M Q , averaged over the 
foreground galaxy population. This is an order-of- 
magnitude higher than the total Ca II mass in the 
Milky Way disk, suggesting that more than 90% of 
Ca II in the Universe is in the circum- and inter- 
galactic environment. 

• The amount of Ca II in the halo is larger for galax- 
ies with higher stellar mass. We find that the mean 
Ca II absorption at fixed virial radius-normalized 
impact parameter is independent of stellar mass. 
This suggests that the Ca II mass in the halo is 
more closely related to the dark matter halo mass 
rather than the total stellar content of galaxies. 

• We find, at fixed stellar mass, more Ca II in the 
halo of star-forming galaxies than quiescent ones 
by a factor of a few. 

• For edge-on galaxies we find that Ca II is preferen- 
tially distributed along the minor axis (in the polar 
direction) . 
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These results provide new constraints on the gas dis- 
tribution around galaxies and bring new elements to our 
understanding of the cosmic baryon cycle. They also 
provide quantitative challenges for theoretical models for 
galaxy formation and evolution: 

o The SFR dependence of the amount of Ca II is simi- 
lar to that of the highly- ionized O VI recently found 
by t he HST COS-halos team ()Tumlinson et al.l 
1201 ID . These results imply star formation, or a 
process correlated with star formation, is respon- 
sible either for injecting Ca II and O VI into the 
halo or keeping a large amount of calcium and oxy- 
gen in the halo in certain physical condition. The 
azimuthal a ngle de pendenc e of Ca II and Mg II 
absorption flBordoloi et al.l I201U fkacprzak et al.l 
120111 : iBouche et al.ll2012tlBordoloi et al.ll2012t Ru- 
bin et al. in prep), which shows that low-ionization 
gas is more abundant and more kinematically vi- 
olent in the polar direction of edge-on galaxies, 
suggests bipolar outflows driven by star forma- 
tion must have played a significant role in inject- 
ing metals into the halo. This is also in agree- 
ment with the ubiquitous outflows found in post- 
starburst /star- forming galaxy spectra at various 
redshifts (e.g., iTremonti et all 120071; I Weiner et al.1 
|200J lkubin et al.ll2010L see lSteidel et al.ll2oTo[ for 
more evidence and a review). Following a period 
of star formation, the fate of Ca II ions is unclear. 

In contrast to the strong SFR dependence of metal 
lines, there appears to be no statistically significant 
difference between the Lya absorption by neutral 
hydrogen around star -forming galaxies and that 
around quiescent ones ()Thom et al.ll2012() , suggest- 
ing neutral hydrogen and metals in the halo could 
have very different origins or be distributed sepa- 
rately in phase space. 

o The amount of Ca II in the halo is correlated with 
stellar mass. This dependence is shallower than lin- 
ear, such that the total Ca II mass is more closely 
related to dark matter mass than stellar mass. This 
may have the same cause as the mass-metallicity re- 
lation - the metallicity of interstellar gas decreases 



with stellar mass because metals are more easily 
expelled from low-mass galaxies due to their shal- 
lower gravitational potential. This will result in 
higher metal (halo)-to-stellar mass ratio for lower- 
mass galaxies. It is interesting to note that the 
total O VI mass in the halo shows even less d epen- 
dence on stellar mass (jTumlinson et al.ll2011|) . 

Large and homogeneous surveys of the sky have 
allowed us to probe the distribution of matter in 
low-density environments. From the Sloan Digital 
Sky Survey only, we now have measurements of the 
galaxy-gala xy correlation functi o n from clustering anal- 
yses (e.g., iMasiedi et all 120061 : IjTang et al.l l2012h . the 
galaxy-mass correlation f rom gravitational lensing (e.g., 
Mandelba um et ail I2006D . the galaxy-dust corre lation 
from reddening measurements (|Menard et al.ll20Tol) and, 
from the work presented in this paper, the galaxy-gas 
correlation function obtained by measuring statistical ab- 
sorption by metals. These new observational probes are 
providing us with a more complete description of the 
matter distribution around galaxies. 

This analysis shows the potential of using the ever- 
growing data from large surveys for absorption line stud- 
ies. The automatic continuum estimation method pre- 
sented in lZhu fc Me nard (2012|) and the cross-correlation 
technique presented in this paper are generic and read- 
ily applicable to any large data set from on going and fu- 
ture surveys such as BO SS (Schl eeel et aHl2009l). e BOSS 
(IComparat et al.ll2012D. B igBOSS (ISchlegel et al.ll20"H , 
and PFS (lEllis et al.H20i2T) . 
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